Model‐assisted identification of metabolic engineering strategies for Jatropha curcas lipid pathways

Efficient approaches to increase plant lipid production are necessary to meet current industrial demands for this important resource. While Jatropha curcas cell culture can be used for in vitro lipid production, scaling up the system for industrial applications requires an understanding of how growth conditions affect lipid metabolism and yield. Here we present a bottom‐up metabolic reconstruction of J. curcas supported with labeling experiments and biomass characterization under three growth conditions. We show that the metabolic model can accurately predict growth and distribution of fluxes in cell cultures and use these findings to pinpoint energy expenditures that affect lipid biosynthesis and metabolism. In addition, by using constraint‐based modeling approaches we identify network reactions whose joint manipulation optimizes lipid production. The proposed model and computational analyses provide a stepping stone for future rational optimization of other agronomically relevant traits in J. curcas.     

Expansion of MIR482/2118 by a class‐II transposable element in cotton

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide‐binding site leucine‐rich repeat (NBS‐LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G. barbadense, G. hirsutum and G. raimondii, but not in G. arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target‐gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class‐II DNA transposable element (TE) in G. raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G. barbadense, G. hirsutum and G. raimondii. The ‘GosTE’ involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3‐bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G. barbadense and G. hirsutum, suggesting a dynamic co‐evolution of miR482/2118 and its NBS‐LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.     

UDP‐Api/UDP‐Xyl synthases affect plant development by controlling the content of UDP‐Api to regulate the RG‐II‐borate complex

Rhamnogalacturonan‐II (RG‐II) is structurally the most complex glycan in higher plants, containing 13 different sugars and 21 distinct glycosidic linkages. Two monomeric RG‐II molecules can form an RG‐II‐borate diester dimer through the two apiosyl (Api) residues of side chain A to regulate cross‐linking of pectin in the cell wall. But the relationship of Api biosynthesis and RG‐II dimer is still unclear. In this study we investigated the two homologous UDP‐D‐apiose/UDP‐D‐xylose synthases (AXSs) in Arabidopsis thaliana that synthesize UDP‐D‐apiose (UDP‐Api). Both AXSs are ubiquitously expressed, while AXS2 has higher overall expression than AXS1 in the tissues analyzed. The homozygous axs double mutant is lethal, while heterozygous axs1/+ axs2 and axs1 axs2/+ mutants display intermediate phenotypes. The axs1/+ axs2 mutant plants are unable to set seed and die. By contrast, the axs1 axs2/+ mutant plants exhibit loss of shoot and root apical dominance. UDP‐Api content in axs1 axs2/+ mutants is decreased by 83%. The cell wall of axs1 axs2/+ mutant plants is thicker and contains less RG‐II‐borate complex than wild‐type Col‐0 plants. Taken together, these results provide direct evidence of the importance of AXSs for UDP‐Api and RG‐II‐borate complex formation in plant growth and development.     

Leucine‐rich repeat receptor‐like kinase II phylogenetics reveals five main clades throughout the plant kingdom

Receptor‐like kinases (RLKs) represent the largest group of cell surface receptors in plants. The monophyletic leucine‐rich repeat (LRR)‐RLK subfamily II is considered to contain the somatic embryogenesis receptor kinases (SERKs) and NSP‐interacting kinases known to be involved in developmental processes and cellular immunity in plants. There are only a few published studies on the phylogenetics of LRR‐RLKII; unfortunately these suffer from poor taxon/gene sampling. Hence, it is not clear how many and what main clades this family contains, let alone what structure–function relationships exist. We used 1342 protein sequences annotated as ‘SERK’ and ‘SERK‐like’ plus related sequences in order to estimate phylogeny within the LRR‐RLKII clade, using the nematode protein kinase Pelle as an outgroup. We reconstruct five main clades (LRR‐RLKII 1–5), in each of which the main pattern of land plant relationships re‐occurs, confirming previous hypotheses that duplication events happened in this gene subfamily prior to divergence among land plant lineages. We show that domain structures and intron–exon boundaries within the five clades are well conserved in evolution. Furthermore, phylogenetic patterns based on the separate LRR and kinase parts of LRR‐RLKs are incongruent: whereas the LRR part supports a LRR‐RLKII 2/3 sister group relationship, the kinase part supports clades 1/2. We infer that the kinase part includes few ‘radical’ amino acid changes compared with the LRR part. Finally, our results confirm that amino acids involved in each LRR‐RLKII–receptor complex interaction are located at N‐capping residues, and that the short amino acid motifs of this interaction domain are highly conserved throughout evolution within the five LRR‐RLKII clades.     

Satellite DNA landscapes after allotetraploidization of quinoa (Chenopodium quinoa) reveal unique A and B subgenomes

If two related plant species hybridize, their genomes may be combined and duplicated within a single nucleus, thereby forming an allotetraploid. How the emerging plant balances two co‐evolved genomes is still a matter of ongoing research. Here, we focus on satellite DNA (satDNA), the fastest turn‐over sequence class in eukaryotes, aiming to trace its emergence, amplification, and loss during plant speciation and allopolyploidization. As a model, we used Chenopodium quinoa Willd. (quinoa), an allopolyploid crop with 2n = 4x = 36 chromosomes. Quinoa originated by hybridization of an unknown female American Chenopodium diploid (AA genome) with an unknown male Old World diploid species (BB genome), dating back 3.3–6.3 million years. Applying short read clustering to quinoa (AABB), C. pallidicaule (AA), and C. suecicum (BB) whole genome shotgun sequences, we classified their repetitive fractions, and identified and characterized seven satDNA families, together with the 5S rDNA model repeat. We show unequal satDNA amplification (two families) and exclusive occurrence (four families) in the AA and BB diploids by read mapping as well as Southern, genomic, and fluorescent in situ hybridization. Whereas the satDNA distributions support C. suecicum as possible parental species, we were able to exclude C. pallidicaule as progenitor due to unique repeat profiles. Using quinoa long reads and scaffolds, we detected only limited evidence of intergenomic homogenization of satDNA after allopolyploidization, but were able to exclude dispersal of 5S rRNA genes between subgenomes. Our results exemplify the complex route of tandem repeat evolution through Chenopodium speciation and allopolyploidization, and may provide sequence targets for the identification of quinoa's progenitors.     

Ptr1 evolved convergently with RPS2 and Mr5 to mediate recognition of AvrRpt2 in diverse solanaceous species

The Ptr1 (Pseudomonas tomato race 1) locus in Solanum lycopersicoides confers resistance to strains of Pseudomonas syringae pv. tomato expressing AvrRpt2 and Ralstonia pseudosolanacearum expressing RipBN. Here we describe the identification and phylogenetic analysis of the Ptr1 gene. A single recombinant among 585 F2 plants segregating for the Ptr1 locus was discovered that narrowed the Ptr1 candidates to eight nucleotide‐binding leucine‐rich repeat protein (NLR)‐encoding genes. From analysis of the gene models in the S. lycopersicoides genome sequence and RNA‐Seq data, two of the eight genes emerged as the strongest candidates for Ptr1. One of these two candidates was found to encode Ptr1 based on its ability to mediate recognition of AvrRpt2 and RipBN when it was transiently expressed with these effectors in leaves of Nicotiana glutinosa. The ortholog of Ptr1 in tomato and in Solanum pennellii is a pseudogene. However, a functional Ptr1 ortholog exists in Nicotiana benthamiana and potato, and both mediate recognition of AvrRpt2 and RipBN. In apple and Arabidopsis, recognition of AvrRpt2 is mediated by the Mr5 and RPS2 proteins, respectively. Phylogenetic analysis places Ptr1 in a distinct clade compared with Mr5 and RPS2, and it therefore appears to have arisen by convergent evolution for recognition of AvrRpt2.